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Abstract 


The three-nucleon bound state problem is studied employing a nucleon- 
nucleon potential obtained from a basic quark-quark interaction in a five- 
channel Faddeev calculation. The obtained triton binding energy is compa¬ 
rable to those predicted by conventional models of the NN force. 
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I. INTRODUCTION 


During the last decade the development of quark-model based interactions for the 
hadronic force has led to nucleon-nucleon (NN) potentials that provide a fairly reliable 
description of the on-shell data. Several models including quark-degrees of freedom have 
been used to study the NN interaction [|I| and also the baryon spectra p[. Among them, 
the chiral quark cluster model is the only one that pursued a simultaneous understanding of 
different low-energy phenomena based on a unique quark-quark interaction. This model is 
able to provide a quantitative description of the NN scattering Q and bound state problems 


Nevertheless, quark-model based NN interactions have not been often used to study few- 
body systems. There might be two different reasons for that. First of all, most of those inter¬ 
action models for the two-nucleon system needed to be supplemented with meson-exchange 
potentials between the baryons to obtain a reasonable description of the experimental data 
loosing in this way their quark-based character. Secondly, other quark-model based in¬ 
teractions were primarily designed to describe the baryon spectra 0] , but lead to unrealistic 
results when they are applied to the two-nucleon system p. 

In this work we want to perform a study of the triton bound state making use of a nonlocal 
NN potential fully derived from quark-quark interactions. The model has been previously 
utilized for investigations of three-body systems {NNN, NNA, A^AA, and AAA), making 
more emphasis on the mass ordering of possible bound states of these systems than on the 
binding energy values 0. In the present work the full nonlocal NN potential will be em¬ 
ployed as it follows from the application of the resonating group method (RGM) formalism. 
This method allows, once the Hilbert space for the six-body problem has been fixed, to treat 
the inter-cluster dynamics in an exact way. Thereby, nonlocalities are generated, reflecting 
the internal structure of the nucleon, which translate into specihc off-shell properties of the 
resulting NN potential. 

Indeed the relevance and/or necessity of considering the nonlocal parts of NN poten¬ 
tials in realistic interactions is still under debate. Over the past few years several studies 
have appeared in the literature which stress the potential importance of nonlocal effects 
for the quantitative understanding of few-body observables and, specifically, for the triton 
binding energy |IU|-p!B[]. However, the majority of these investigations [pH|-jT5| explore only 


nonlocalities arising from the meson-exchange picture of the NN interaction. 

The nonlocalities generated in a quark-model derivation of baryonic potentials may play 
a relevant role for the case of the three-nucleon bound state. It has been argued that the 
assumptions associated with meson-exchange models sharply limit the nature of the off- 
shell properties of those potentials, once the on-shell matrix elements are constrained to £t 
the NN data 
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Therefore, it is very interesting to investigate the off-shell features of 
potentials derived from a quark-model. Some preliminary studies in this direction have been 
done in Ref. [^. However, there only the short-range part of the interaction is obtained 
by means of quark-model techniques. The intermediate and long-range part is described 
by “standard” meson-exchange between baryons. Accordingly, that work allows only very 
limited conclusions with regard to effects of the quark substructure. 

The triton binding energy is obtained from a Faddeev calculation. We restrict ourselves 
to the standard five channel case; i.e. we consider only the and ^Si — ^Di NN partial 
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waves, those which provide the bulk contribution to the three-nucleon binding energy. The 
three-body Faddeev equations will be solved in momentum space making use of separable 
hnite-rank expansions of the two-body interactions. 

The paper is organized as follows. In Section II we introduce the basic quark-quark 
Hamiltonian and we describe the method to obtain the RGM NN interaction. In Section 
III we provide details about the hnite-rank expansions of the quark-model based potentials 
which enter in the Faddeev calculations of the triton binding energy and we present results 
for the three-body system. Finally, some concluding remarks are provided in Sec. IV. 


II. QUARK-MODEL BASED NN POTENTIAL 


In recent years a chiral quark cluster model for the NN interaction has been developed. 
This model has been widely described in the literature therefore we will only 

briehy summarize here its most relevant aspects. It contains, as a consequence of chiral 
symmetry breaking, a pseudoscalar and a scalar exchange between constituent quarks coming 
from the Lagrangian 


Lch = gchF{q'^)^{a + ■ 7f)T , 


( 1 ) 


where F{q^) is a monopole form factor 


F{q^) = 


A2 

X 


A2 -K 


( 2 ) 


Ay determines the scale of chiral symmetry breaking, being bound between 1 GeV and 600 


MeV 1^. The chiral coupling constant Qch is chosen to reproduce the experimental nNN 
coupling constant. 

/^From the above Lagrangian a pseudoscalar {PS) and a scalar (5) potential between 
quarks can be easily derived in the non-relativistic approximation. 


Kf (?) = 


9ch 






A^ -|- g 2 rripg + q^ 




( 3 ) 


A^ 

rS/^ _ „2 X 


Vii {^ = -9ch 


Ay -|- q^ rrig -|- 


( 4 ) 


Using the range of values for Ay given above yields a V — A mass difference due to 
the pseudoscalar interaction between 150 and 200 MeV. The rest of the mass difference, up 
to the experimental value, must have its origin on perturbative processes. In the present 
model, this is taken into account through the one-gluon exchange potential ||20||, 


VS^^{^ = ■ A,) 


TT 


gZ 


TT 


4m^ 


1 + 


+ 


vr 










( 5 ) 


where the A’s are the color Gell-Mann matrices and Ug is the strong coupling constant. 
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For the present study we make use of the nonlocal NN potential derived through a 
Lippmann-Schwinger formulation of the RGM equations in momentum space El. The 
formulation of the RGM for a system of two baryons, Bi and B 2 , needs the wave function 
of the two-baryon system constructed from the one-baryon wave functions. The two-baryon 
wave function can be written as: 


where A is the antisymmetrizer of the six-quark system, x(P) is the relative wave-function 
of the two clusters, is the internal spatial wave function of the baryon Bi, are 

the internal coordinates of the three quarks of baryon Bi. x^B 2 denotes spin-isospin wave 
function of the two-baryon system coupled to spin (S) and isospin (T), and, hnally, ^c[2^] is 
the product of two color singlets. 

The dynamics of the system is governed by the Schrodinger equation 


{ n - ET )\^>=0 => Et )\^ >= 0 , 


(7) 


where 


N ->2 


^ 2m, 


Kj 


( 8 ) 


with Tc.m. being the center of mass kinetic energy, Vij the quark-quark interaction described 
above, and m, the constituent quark mass. 

Assuming the functional form 


4>b(p) 



(9) 


where b is related to the size of the nucleon quark core, Eq. can be written in the 
following way, after the integration of the internal cluster degrees of freedom, 

x(P) + J {Vd{P, P) + Wl{P, Pi)) x{p)dp = 0 (10) 



Vd{P, Pi) is the direct RGM kernel and Wl{P, Pi) is the exchange RGM kernel, composed 
of three different terms 


WLiP, Pi) = Tl{P, P) + VLiP, Pi) + {E + Ein)KL{P, Pi ), ( 11 ) 

where Ein is the internal energy of the two-body system, TB{P,Pi) is the kinetic energy 
exchange kernel, VL{P,Pi) is the potential energy exchange kernel and KL{P,Pi) is the 
exchange norm kernel. Note that if we do not mind how VD{P,Pi) and WiiP^Pi) were 
derived microscopically, Eq. ( pUj) can be regarded as a general single channel equation of 
motion with including energy-dependent non-local potential. Vd{P, Pi), which contains the 
direct RGM potential, and Wl{P', Pi), which contains the exchange RGM potential coming 
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from quark antisymmetry, constitute our energy-dependent nonlocal potential. In our case 
Ein = 2mN what makes our potential almost energy-independent, because the center of 
mass energy of the two-body system E is much smaller than the internal energy Ei^- 

The potential yields a fairly good reproduction of the experimental data up to laboratory 
energies of 250 MeV. For a correct description of the phase shift it is necessary to take 
into account the coupling to the NA channel 
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that provides an isospin-dependent 
mechanism generating the additional attraction in this channel. This is implemented in our 
calculation generalizing Eq. ([T0|) to a coupled channel scheme. It implies a modihcation of 
Eq. o with an additional term which contains the NN —>■ NA coupling. The parameters 
used are summarized in Table |. In Table y we present the results for the low-energy 
scattering data and the deuteron properties of the present model together with values of 
some standard NN potentials m |2^ and experimental data. It is known that a charge 
symmetry breaking term should be included in the interaction if one wants to reproduce 
those quantities simultaneously |^. This is taken into account by a slight modification 
of the chiral coupling constant to reproduce the deuteron and the low energy scattering 
parameters (see Table |^. We also show, in Figs. Qg and 1^, the and ^Si — ^Di phase 
shifts and the mixing parameter ei in comparison to results from phase-shift analyses 


III. TRITON BINDING ENERGY 


The triton binding energy is obtained by means of a Faddeev calculation using the NN 
interaction described above. We perform a so-called hve-channel calculation, i.e., we use 
only the ^Sq and — ^Di NN partial waves as input. Note that since in our model 
there is a coupling to the NA system, as explained above, a fully consistent calculation 
would require the inclusion of two more three-body channels. However, their contribution 
to the 3N binding energy is known to be rather small and therefore we neglect them 
for simplicity reasons. 

To solve the three-body Faddeev equations in momentum space we first perform a sepa¬ 
rable finite-rank expansion of the NN{—NA) sector utilizing the EST method Such a 
technique has been extensively studied by one of the authors (J.H.) for various realistic NN 
potentials and specihcally for a model that also includes a coupling to the NA system 
|30| . In those works it was shown that, with a separable expansion of sufficiently high rank, 
reliable and accurate results on the three-body level can be achieved. In the present case 
it turned out that separable representations of rank 6 - for — (^Dq) and for ^Si — ^Di 
- are sufficient to get converged results. The set of energies used for the EST separable 
representations are listed in Table B We refer the reader to Refs. 


3T| for technical 


details on the expansion method. The quality of the separable expansion on the NN sector 
can be seen in Fig. |l|, where we show phase shifts for the original nonlocal potential and for 
the corresponding separable expansion. Evidently, the phases are almost indistinguishable. 

Results for the triton are summarized in Table IV. It is reassuring to see that the predicted 
triton binding energy is comparable to those obtained from conventional NN potentials, such 
as the Bonn or Nijmegen models. Thus, our calculations show that quark model based NN 
interactions are dehnitely able to provide a realistic description of the triton. The results 
also give support to the use of such interaction model for further few-body calculations. 
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One should not forget at this point that the number of free parameters is greatly reduced 
in quark-model based NN interactions like the present one. In addition, the parameters 
are strongly correlated by the requirement to obtain a reasonable description of the baryon 
spectrum. 


IV. CONCLUSIONS 


We have calculated the three-nucleon bound state problem utilizing a nonlocal NN 
potential derived from a basic quark-quark interaction. This potential was generated by 
means of the resonating group method so that nonlocalities resulting from the internal 
structure of the nucleon were preserved. The resulting triton binding energy is comparable 
to those obtained from conventional NN potentials. 

In the calculation of the three-nucleon binding energy we have followed the traditional 
approach, namely solving the Faddeev equations with nucleon degrees of freedom. Let us re¬ 
mark, however, that in a more fundamental approach one would impose consistency between 
the treatment of the two- and three-nucleon systems in terms of quark degrees of freedom. 
That, of course, would require the derivation and solution of corresponding three-nucleon 
RGM equations. In such a framework the quark structure of nucleons generates (besides the 
consecutive two-nucleon interactions that are summed up by the Faddeev equations) also 
genuine three-body forces. These forces are of short-ranged nature and they could be signih- 
cant for short-distance phenomena like the high-momentum-transfer part of the charge form 
factor of ^He. Indeed, there have been attempts to explore the effects of such three-body 
forces on the triton binding energy. In a simple model based on a single one-gluon exchange 
3^ the three-body exchange kernels have been evaluated. An estimation provided in this 


reference suggests that those three-nucleon forces could yield additional binding in the order 
of 0.2 MeV. If this is the case then those effects would be still small enough to guarantee 
that the approach we followed in our study is sufficiently accurate for an exploratory cal¬ 
culation. However, one has to keep in mind that the estimation in Ref. was done only 
in perturbation theory, and by means of a zeroth order three-nucleon wave function with a 
series of htted parameters. Thus, for the future, a more rehned and consistent treatment of 
the three-nucleon problem within the quark picture is certainly desirable in order to allow 
for reliable conclusions on this issue. 
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TABLES 

TABLE 1. Quark model parameters. The values in brackets are used for a correct description 
of the deuteron. 


mq(MeV) 

313 

6^(fm) 

0.518 

O-s 

0.4977 

9ch 

6.60(6.86) 

m5(fm“^) 

3.400 

mp5(fm“^) 

0.70 

Ax(fm-i) 

4.47 


b is the parameter of the harmonic oscillator wave function used for each quark r]{x) = 

)' 


. Trb' 




TABLE II. NN properties 




Quark model 

Nijm II 

22 

Bonn B | 

21 

Exp. 

Low-energy scattering parameters 


as (fm) 

-23.759 

-23.739 



-23.750 


-23.74±0.02 


Vs (fm) 

2.68 

2.67 



2.71 


2.77±0.05 


at (fm) 

5.461 

5.418 



5.424 


5.419±0.007 


n (fm) 

1.820 

1.753 



1.761 


1.753±0.008 

Deuteron properties 

ed (MeV) 

-2.2242 

-2.2246 



-2.2246 


-2.224575 

Pd (%) 


4.85 

5.64 



4.99 


- 

Qd (fm^ 


0.276 

0.271 



0.278 


0.2859±0.0003 

As (fm“ 

1/2) 

0.891 

0.8845 



0.8860 


0.8846±0.0009 

Ad! As 


0.0257 

0.0252 



0.0264 


0.0256±0.0004 


TABLE III. Expansion (lab) energies (in MeV) used in the EST representations of the 
quark-model potential, refers to the deuteron binding energy. Q is the boundary condition 
chosen for the angular momentum Q of the initial state pO|,pT[. 


Partial wave (E^,Q) 

^So{NN) -^Do{NA) (TO) (50,0) (300,0) (-20,0) (-20,2) (-50,0) 

^Si - ^Di €d (100,0) (175,2) (300,2) (-50,0) (-50,2) 























TABLE IV. Properties of the three-nucleon bound state. 



Quark model 

Nijm II 

Bonn B 

Eb (MeV) 

-7.72 

-7.65 

-8.17 

Ps (%) 

91.49 

90.33 

91.35 

Ps' (%) 

1.430 

1.339 

1.368 

Pp (%) 

0.044 

0.064 

0.049 

Pd (%) 

7.033 

8.267 

7.235 
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FIGURES 


FIG. 1. ^So NN phase shift. The solid line is the result for the nonlocal quark-model potential. 
The dotted line shows the result of the separable representation of the nonlocal quark-model 
potential. The squares, diamonds and triangles are the experimental data taken from Refs. |24|, 
, and [^1, respectively. 


FIG. 2. Phase shifts for the ^Si and ^Di partial waves. Same description as in Fig. |^. 
FIG. 3. Mixing parameter ei. Same description as in Fig. |l|. 
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